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Abstract 

Next to leading order effective field theory calculations are performed for 1 5'o 
NN scattering using subtractive renormalization procedure. One pion exchange 
and contact interaction potentials are iterated using Lippman-Schwinger equa- 
tion. Satisfactory fit to the Nijmegen data is obtained for the momenta up to 
300 MeV in the centre of mass frame. Phase shifts are also compared with the 
results of KSW approach where pions are included perturbatively. 
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I. INTRODUCTION 



There has been much recent interest in applications of the chiral perturbation theory ap- 
proach to processes involving an arbitrary number of nucleons introduced by Weinberg [l].^] • 

For processes involving more than one nucleon Weinberg's power counting is applied to 
the potential rather than to the scattering amplitude. For n-nucleon processes the potential 
is defined as a sum of n-nucleon irreducible time-ordered perturbation theory diagrams. The 
amplitude is obtained by solving Lippmann-Schwinger equation (or Schrodinger equation). 

Iteration of the potential via the Lippmann-Schwinger equation leads to divergences. 
One could try to regularize the potential and include counter-terms, but due to the non- 
renormalizability of the theory one would have to include an infinite number of (counter-)terms 
with more and more derivatives. Hence, one has either to exactly solve the equation and after 
subtract divergences explicitly, or otherwise draw all relevant diagrams, subtract them and sum 
up. There are not any equations for subtractively renormalised amplitudes. Fortunately there 
exist practically more powerful approach, cut-off EFT j4|. Divergences appearing in diagrams 
can be regulated using (sharp or smooth) cut-off regularization. One can keep cut-off parameter 
of the order of the mass of lightest integrated particle and fit coupling constants to the exper- 
imental data as was done in ref. ||. Cut-off EFT is equivalent to subtractively renormalised 
EFT up to the order of considered accuracy 

Technically more convenient approach based on a new systematic power counting has been 
suggested by Kaplan, Savage and Wise (KSW counting) in || . Pions are included perturbatively 
within this new scheme. A similar power counting with perturbative pions has been suggested 



by Lutz ||. As was mentioned in |T0[ the expansion parameter in KSW counting is rather 
large. It raises a question of the usefulness of suggested expansion. Arguments supporting as 
well as criticising the perturbative inclusion of pions have been given in the literature (see [|TT] 



and citations therein. ]TT] also contains references to lots of papers devoted to EFT study of 
NN interaction problem). 

In Weinberg's power counting the one pion exchange potential is of leading order and hence 
it has to be iterated via the Lippmann-Schwinger equation. KSW counting suggests that 
contributions containing multiple iterations of the one-pion exchange potential are higher order 
and hence need not be included in low order calculations. 

In this letter the results of the application of subtractive renormalization procedure to the 
Weinberg's approach to 1 5'o NN scattering problem are briefly discussed (the details will be 
given in separate paper). Leading order (contact interaction plus one-pion exchange) plus next- 
to-leading order contact interaction potential is iterated using Lippman-Schwinger equation. 

According Weinberg's power counting the leading order potential for 1 Sq NN scattering is 
the following: 



Vo(p,P , ) = C' + K(p,P , ) (1) 



where 



2/2 (q 2 + ml) 8vr/2 

C is determined by contact interaction terms in the effective Lagrangian ||, q = p — p', 
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FIG. 1. Diagrams defining A^, Xp and Ge- Solid and dashed lines correspond to nucleon and 
pion respectively. Small circle in the diagram defining Xp indicates that factor C is excluded from that 
vertex. 

= 1.25 is the axial coupling constant, m n = 140 MeV is the pion mass and f^ = 132 MeV 
is the pion decay constant. 

Substituting Vq into the Lippman-Schwinger equation one finds for the Feynman amplitude 
13: 



C[Xpf 
L - CGi 



(3) 



where the quantities A n , Xp an d Ge correspond to diagrams drawn in FIG. 1. A n and Xp 
are finite. The first diagram contributing to Ge is linearly and the second logarithmically 
divergent, all other diagrams are finite. As A contains divergences it is necessary to regularize 
and renormalize it. 

Let us introduce a cut-off into the potential: 



(p 2 + A 2 ) 2 (p' 2 + A 2 



+ K(p,p) 



(4) 



where A is a cut-off parameter. 



Substituting into the Lippmann-Schwinger equation one gets the following expression 



for the amplitude: 



xA — ^A 7 



C 



Xr 



l-CG E 



(5) 



where A^ is unchanged, Xp an d G E are given by the same diagrams in FIG.l with C replaced 
by ^ - 



Note that Xp = X P 



A _ „ + o(l/A) and G E = a x A + a 2 lnA + G f E + 0(1/A), where G f E is finite 
and A- independent. It is straightforward to find a\ and a 2 by calculating analytically first two 
diagrams contributing to G E . 
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FIG. 2. [xp(0)] plotted versus momentum in MeV. The solid and dashed lines correspond to the 
real and imaginary parts respectively. 



One can numerically calculate A n and A A for different values of A. Using these values of 
A„ and A A one finds: 



1 



A — Ak 



B 1 (p)A + B 2 (p) InA + B 3 (p). 



(6) 



Although B 3 (p) actually depends on A, for very large A it practically does not and Bi(p) and 
B 2 {p) can be calculated very accurately. Having calculated ai, a 2 , Bi(p) and B 2 (p) one easily 
finds: 



Bi(p) 



(7) 



Xp 



«2 



S 2 (p) 



(8) 



Good agreement between numerical results obtained from (|7|) and (§) suggests that numerical 
calculations are enough accurate. [Xpf obtained this way is plotted in FIG. 2. The same 
quantity was calculated in |12]| using a different method. The agreement between FIG. 2 and 



the results of |I2[ is quite satisfactory. 



Having calculated 



A 

Ap 



(« [xp] 2 ) 



it is not difficult to calculate the sum 



of all 



finite regularised (unregularised) diagrams contributing to [Ge)- 

Phase shifts corresponding to subtractive renormalization procedure can be calculated per- 
forming following steps: 

1. Calculate A n , calculate A A for very large A (A ~ 10 5 MeV has been taken); 

2. Calculate 
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FIG. 3. Phase shifts in degrees versus centre of mass momentum. Solid line corresponds to Ni- 
jmegen phase shift data. Long-dashed line is the leading order EFT result. Dotted and short-dashed 
lines are the results of next-to-leading order calculations. The dotted and small-dashed lines correspond 
to non-perturbative and non-perturbative inclusion of next-to-leading order potential C2 (p 2 + p'j /2. 




W(p) = -- d + G 



A^ A 7 



i A (p) - lib) 



(9) 



where I A (p) and l£(p) correspond to the first and second divergent diagrams contributing to 
respectively. 

3. Calculate two divergent diagrams contributing to analytically, subtract the second 
logarithmically divergent diagram at p 2 = —fi 2 and the first linearly divergent diagram at p 2 = 
—v 2 [= — /t/i 2 ) (one is free to subtract these two diagrams at different points). In calculations of 
this work fx = 240 MeV and v = 800 MeV were taken (k = 100/9). Note that the large value 
of the second subtraction point (800 MeV) does not affect the power counting. 

Subtracted diagrams give: 



g\ fm n M\ 



2f 2 V 4tt J 



In 1 



mi 



2/i 



m, 



and 



jR = M(u + ip) 



4. Calculate renormalized amplitude: 

A R = At, 



4:71 



(10) 



:n) 



W{p) + I R (p) + I*{p) 



12) 
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5. Fit C so as to reproduce phase shifts for very low energies (p ~ 5 MeV). One gets 
C = -1/(17348 MeV 2 ). 

6. Calculate phase shifts. Comparing phase shifts calculated for A = 1 x 10 5 MeV, 2 x 
10 5 MeV, 3 x 10 5 MeV, 4 x 10 5 MeV one sees that A R is practically cut-off independent for very 
large A. 

These leading order phase shifts are plotted in FIG. 3. 

As was argued in |T2[ Ci coupling constant of the part of the next-to-leading order potential 

V 1 = C 2 (p 2 + p' 2 ) /2 (13) 

is larger than assumed by Weinberg's power counting. This suggests that contributions of this 
part of the next-to-leading order potential are larger than of others. This observation simplifies 
the actual next-to-leading order calculations significantly. It is possible to sum up all diagrams 
obtained by iterating Lippman-Schwinger equation for the potential V = Vq + Vi, where Vo 
and V\ are given by (H) and (0) respectively. The result obtained in |12] using dimensional 
regularization reads: 

A 1 (p) = A 1t -- ^£ — 1 (14) 



C - a A MC 2 + C 2 p 2 



G } 



As x p and Ge are already calculated (although these quantities were calculated using cut-off 
regularization, as was mentioned above the effect of finite cut-off is extremally small and hence 
negligible), it is straightforward to apply the same subtraction scheme to ([14]), fix C 2 so as to 
fit the data for p = 10 MeV and calculate the phase shifts. The results are plotted in FIG. 3. 



One can see that agreement with Nijmegen phase shift data |13] is quite good. Obtained value 
C*2 = 1/(8.5 x 10 9 MeV 4 ) is indeed larger than assumed by Weinberg's counting. 

One also can include V\ perturbatively expanding (|14|) in C 2 and retaining only first two 
terms. Calculated phase shifts are also plotted in FIG. 3. Comparing the phase shifts for 
perturbative and non-perturbative inclusion of V\ one sees that deviation is small for momenta 
up to 300 MeV. This suggests that higher order contributions of V\ are small as is expected 
from power counting arguments. 

On the other hand it is straightforward to include pions perturbatively (as was done in 
TJJ) following ideas of KSW counting. Using Those values for C and C 2 which were obtained 



above one gets negative phase shifts, this suggests that the higher order contributions of the 
one pion exchange potential are by no means small. Note that subtraction scheme used here 
is different from the one considered in ||14|| . In [14| finite diagrams were subtracted as well. 



Repeating the calculations of this work with non-perturbatively included one pion exchange 
potential using PDS scheme and substituting the "best fit" parameters from [§] one gets again 
negative phase shifts. This comparison of the phase shifts calculated using perturbative and 



non-perturbative inclusion of pions supports the conclusion of that the consistency of the 



perturbative inclusion of pions in EFT approach to the A^A^-scattering problem is questionable. 
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